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Soft Physics in STAR 

G. Van Buren a for the STAR Collaboration* 

a Brookhaven National Lab, Upton, New York 11973-5000, USA 

The STAR Experiment at RHIC is well-suited to making measurements of particle yields 
from relativistic nuclear collisions in the low transverse momentum, or "soft" , regime. We 
present preliminary results on measurements of ir°, tt~ , ir + , K~, K + , K®, p, K*°+K*°, fo, 
p, p, 0, A, A, E~, S + , n~, and fi + from the first two years of physics running at RHIC. 
An abundance of physics topics can be addressed by these measurements, some of which 
are discussed here along with the results. 

1. Introduction 

"The principles of thermodynamics are devised for giving a phenomenological account of 
the gross behavior of macroscopic physical systems, ...[interpreting] the average behavior 
of systems of many degrees of freedom." [] The ongoing search for a system of nuclear 
matter which can be described by a thermodynamic quark gluon plasma phase (QGP) 
must therefore involve a study of nuclear matter in bulk. For this reason, heavy ions 
with significant amounts of nuclear matter are used in laboratory experiments to find 
and study QGP. This also implies that QGP formation affects some bulk portion of the 
outgoing matter, most of which is in the "soft", low-pr region of phase space. 

The STAR Experiment at RHIC has measured a rich variety of particle species and 
portions of their soft spectra in high energy Au+Au and p+p collisions [ D], 0, [J], £|, [5], [(], |7|, 
§. Because production, annihilation, and rescattering processes differ among the particle 
species, their abundances and spectra reflect on different aspects of the collision evolution. 
We present here a discussion of these aspects as elucidated by our measurements. We will 
follow an outline of the physics topics along the evolution from initial collisions to final 
reinteractions of the products. 



2. Inclusive Multiplicity and Mean Transverse Momenta 

Improved statistics in the second year of physics running at RHIC have given STAR the 
ability to measure unidentified charged hadron transverse momentum (jpx) spectra from 
0.2 GeV/c to ~12 GeV/c. These spectra are shown in Fig. [I] as a function of centrality for 
Au+Au collisions. They have been corrected for decays and backgrounds through detailed 
simulation [ [|. The corrections are px and multiplicity dependent, but nowhere exceed 
20%. From these inclusive spectra we can measure total yields as well as (pr)- This is done 



*For the full author list and acknowledgements, see Appendix "Collaborations" of this volume. 
^R. Tolman, The Principles of Statistical Mechanics, Oxford University Press (1938). 
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by fitting the spectra with a power law of the form C(l+pr/po) n [ 0- The pr-integrated 
multiplicity requires extrapolation of the power law at the low-p^ end, which comprises up 
to 21% of the yields. For the 5% 
most central events, our preliminary re- 
sults indicate that the average charged 
hadronic multiplicity (dN c h/drj) increases 
from 563±39 to 687±69 between ^s NN 
= 130 and 200 GeV. (p T ) is then ob- 
tained from the same power law, giving 
0.520±2.3% GeV/c at = 130 GeV 

and 0.524±2.3% GeV/c at 200 GeV. For 
both (dN ch /drj) and (px) measurements 
the errors are dominated by systematics 
mostly common between the two energies 
and do not reflect possible relative varia- 
tions between each other. 

To put these measurements in context, 
a parameterization of the Color Glass 
Condensate (CGC) model 
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fit to central A+A (p T ) at -ys NN 
130 GeV in Fig. | [ |T|. While the CGC 



Figure 1: Transverse momentum spectra of hadrons 
from y/s NN = 200 GeV Au+Au over several central- 
izes. Only statistical errors are shown. See Ref. [ | 
for systematic errors and more details. 

parameterization can fit these two data 

points, the 200 GeV central Au+Au (pt) datum falls well below the fit from the lower 
energies, despite anticipation that the model is most applicable in high energy central 
A+A collisions. The CGC model has had numerous other successes 
supplemental considerations to congrue with these data [ [TJ . 

Also evident in Fig. Q are the differences 
in collision energy systematics for (p?) be- 
tween A+A, e + +e~ and p(p)+p data. This 
implies that A+A is not just a simple su- 
perposition of nucleon+nucelon collisions [ 



13] , but may require 



0.6r 



'0.55- 



0.5- 



e +e 
□ OPAL 

A AMY 
O TASSO 
* Mark 



A+A (central) 

• NA4S 

★ STAR Preliminary 



14 



and argues against a universality in 
particle production [[HJ suggested from to- 
tal particle yields [ [U 



3. Baryon Number 

Conserved quantities must survive any 
evolution. The conserved net-baryon num- 
ber is established from the incoming nu- 
clei, but becomes spread over the rapidity 
interval by initial interactions (often called 
stopping, though not strictly stopped), and 
then additionally smeared by the rescatter- 
ing processes; together these make up the 
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Figure 2: Mean transverse momenta of hadrons 
for central heavy ion, e + +e~, and p(p)+p collisions 
versus energy [ [l2| . Curves are CGC model for A+A , 
phenomenological parameterization for p(p)+p, and 
JETSET calculations (versus thrust axis) for e + +e~. 
STAR data errors are systematics-dominated. 
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transport of baryon number from the nuclei to their final rapidities. The midrapidity 
baryon and antibaryon yields thus have folded into them these transport distributions 
along with components from baryon-antibaryon pair production (and annihilation, which 
may also contribute to transport) processes. 

A characterization of the baryon-anti- 
baryon differences is made in Fig. |3|. This 
plot shows the antibaryon to baryon ra- 
tios at midrapidity for various species in 
central heavy ion collisions. It is clear 
from the plot that the ratios tend to- 
wards one with increasing strangeness con- 
tent, and with increasing energy. The 
ratios are not corrected for feeddown, 
and the p and p are essentially inclu- 
sive. These ratios indicate that the yield 
of baryons resulting from pair production 
processes has begun to dominate over those 
from transport at RHIC energies: the ra- 
tio of those contributions from inclusive 
(anti)protons in central Au+Au is ~2.5 
at ^s NN = 130 GeV, derived from inclu- 
sive p/p=0.72±0.05(syst.), and ~3.5 at 
200 GeV, from p / p=0. 7 '8±0. 05 (syst.). In- 
clusion of strange quarks in the baryons ap- 
pears to dilute the contributions of transport, enough to make the Q /Q~~ ratio consistent 
with 1 at RHIC; this is not the case at SPS energies where the contributions from transport 
are still evident in the Q yields despite an initial net strange quark number of zero. 

STAR has also measured most of these antibaryon to baryon ratios as a function of 
centrality and sees only a slight rise in the antiproton to proton ratio for peripheral 
collisions, and no significant changes for the other ratios [ |3], [l7j • This implies little 
change in the relevance of pair production processes for baryon and antibaryon production 
despite the changing collision geometry and participants. 

A definitive measure of the baryon number transported to midrapidity comes from the 
difference between yields of baryons and antibaryons at midrapidity. This is represented 
by the net inclusive proton number shown in Fig. §(left). This number is clearly decreasing 
with increasing collision energy for large heavy ion collisions, while the yield of produced 
particles at midrapidity (represented in Fig. U(left) by the negative hadron yield) is mono- 
tonically increasing. This demonstrates the decreasing significance of baryon transport 
contributions to midrapidity yields and is not unexpected as the rapidity interval from 
the original nuclei is stretched. It is interesting that the baryon number transport is in 
fact still nonzero despite the interval of ~5 units in rapidity at 130 GeV (also true for ~6 
rapidity units at 200 GeV). Some understanding of this may come from the asymmetric 
q versus q distributions even at low momentum fraction x=p q q/p^ within the incoming 
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Figure 3: Midrapidity antibaryon to baryon ratios 
are shown for various species and at several energies 
of central heavy ion collisions. Only statistical errors 
are shown. Systematic errors range from ~7% for 
p/p to -15% for in RHIC data. 
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Several model predictions are also shown for the centrality plot, with 



HIJING matching peripheral data best, and HIJING/BB (which is intended to enhance 
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Figure 4. The net protons at midrapidity are shown versus collision energy for central collisions along 
with the negative hadrons (left, statistical errors only), and versus centrality for ^/s NN = 130 GeV along 
with model predictions (right, errors systematics-dominated) . 

the stopping [ |H| ) more closely approaching the central data. RQMD overpredicts the 
transport contributions. 

4. Particle Production 

Long-standing arguments exist for extraordinary hadron yields should a QGP form, 
such as strangeness phase space saturation (enhanced yields over a solely hadronic evo- 
lution) [ |20| and J/tp suppression [ |2lf . While studies of hadron production in the low 
Pt region of phase space where a QGP is likely to manifest are appropriate for these 
searches, it should be understood that experiments often measure yields after any QGP 
has condensed, at the end of all hadrochemical reactions which can potentially hide earlier 
irregularities (equilibration of a hadron gas, for example). 

It is all the more important that we study the particle production to learn what we 
can about the system and its evolution. STAR is well-suited for this task, capable of 
identifying and measuring spectra and yields for a plethora of particle species. Results on 
the following species have been presented at this conference: 7r°, 7r~, 7r + , K~ , K + , K®, 
p, f , p, p, 0, A, A, S", S + , n~, and [ g g, |, §, §, % §. The impressive 

preliminary spectra and fits for several of these species from the y/s NN = 200 GeV data 
can be seen in Fig. |5], and in Table |I] where the rriT reach of the measurements and results 
of the fits are presented. Our observations of the K*°, p, and fo are the first in heavy ion 
collisions. 

One can attempt to obtain a global characterization of particle production by compar- 
ing the yields to a model which describes the bulk properties of the system through a 
limited set of parameters. A thermal, chemical equilibrium model is one such statistical 
model, and we can examine how this model describes our y/s NN = 130 GeV midrapid- 
ity data, mindful of any caveats [ |2"2"| . In the particular model we choose, free param- 
eters include the temperature (T), light quark chemical potential (p g ), strange quark 
chemical potential (p s ), and strangeness saturation factor (7 S , where values less than 
one indicate incomplete saturation) [ j23fl . Shown in Fig. is a fit excluding the K*, 
(ft, 5, and Q results (x 2 /dof = 1.53). From this we can compare the excluded data 
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Figure 5. A variety of preliminary midrapidity spectra from ^/snn = 200 GeV Au+Au and p+p collisions 
in STAR are plotted, including those of tt", K~ , and inclusive p [§, (iT*°+#*°)/2, p, and f [ §, <j> [ §, 
and f2 [ U . Only statistical errors are shown. See respective references for systematic errors and further 
details. All fits are uit exponentials except for Bose-Einstein on tt - and pt Gaussian on p\ the fit results 
can be found in Table 111 for selected centralities. 
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to the outcome of the fit based on the included species, demonstrating how well the 
model matches all of the data. Including all ratios results in % 2 /dof = 1.68. In ei- 
ther fit, the preferred value of 7^ is close to 1. Other chemical equilibrium models have 
also predicted and fit most of these particle yields rather well [ [M], converging on sim- 
ilar values for the thermodynamic parameters (with a temperature range from 165±7 
MeV to 179±9 MeV, for example). While these models may or may not be applica- 
ble, and good fits are neither new (even seen in e + +e~ and p(p)+p data, albeit without 



strangeness saturation) nor absolute proof of equilibrium 
1) the fits indicate approximate 
strangeness phase space satura- 
tion for central collisions, and 2) 
recent lattice QCD calculations 
for non-zero chemical potentials 
have determined a similar temper- 
ature of 172±3 MeV in this re- 
gion of the phase boundary be- 
tween hadronic matter and QGP [ 



25| , our results are provocative: 



-1 

010 



-2 
10 



26 



Table shows the preliminary 
results of spectra measurements 
in STAR from Au+Au at ^s NN 
= 200 GeV for vr", K~, p, 
K*°+K*°, p, f , 0, fi" and U + [ 
|2], |, ||, |j. A preliminary px- 
independent muon contamination 
correction has been applied to the 
pions. More data is necessary to 
tightly constrain a chemical fit. 
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Figure 6: A fit to the STAR midrapidity 15% central 130 GeV 
Au+Au data with a thermal statistical model [ E3I] is shown. 



5. Transverse Momentum Systematics 

Transverse radial flow can be gauged 
through the mass dependence of either inverse 
slope parameters (T) of spectra fits (inversely 
proportional to the steepness of yield fall-off 
in mr) [ p8| , or (p?) [ of identified spec- 
tra. These are shown for various STAR data 
from central Au+Au collisions at 130 GeV in 
Figs. [?| and ||(left). It is important that the 
inverse slope is measured over a large range of 
the produced spectrum to permit valid com- 
parisons, and this holds for the data shown 
in Fig. [7|. Here it can be seen that the rising 
trend versus mass expected from transverse 
radial flow [ ^] is seen for the lighter species, 
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Figure 7: Inverse slope parameters as a function 
of particle mass for SPS [ |7|, || and RHIC are 
shown with statistical error bars, and systematic 
error brackets for RHIC data only. 
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Figure 8. Mean transverse momenta of identified species versus mass (left) for ^/s NN = 130 GeV 
Au+Au collisions along with blast wave model fit (upper curve: Tfo— 100 MeV, (/3t)— 0.55) and zero flow 
curve (lower curve: Tfo— 100 MeV, f3r=0). (pr) of 7r~, K~, andp are mapped against multiplicity (right) 
for y / SNN = 130 GeV Au+Au, and 200 GeV Au+Au and p+p. All errors are systematics-dominated. 



but that the strange baryons do not clearly follow this trend. Further discussion of these 
systematics can be found in Ref. [ |12|j . 

The same behavior is exhibited in Fig. §(left) by the central (pt) values. Fig. ||](right) 
shows how the (pr) develops as a function of centrality for 7r~, K~ , and p at RHIC. Going 
from p+p data, through peripheral to central Au+Au data shows a smooth rise consistent 
with the onset and buildup of transverse radial flow. 

A blast wave [ [30] fit to the central 7r~ , K~ , and p spectra provides a transverse radial 
flow velocity of (/3 r )~0.55 at 130 GeV, and ~0.60 at 200 GeV. These flow velocities 



appear higher than the values of ~0.50 obtained at the SPS [|3~If, indicating stronger flow. 
However, some deviation from the global blast wave fit is seen in the multistrange baryons 
at RHIC [ffl. One additional note is that the fit kinetic freezeout temperature parameters 
show at most a mild decrease going from SPS (T fo ~110 MeV) to RHIC (T fo ~100 MeV) 
energies. 



6. In-Medium Decay 

One exciting prospect is the possibility to probe the time scale of hadronic rescattering 
via in-medium decay. The principle here is that the daughters of a particle which decays 
during the rescattering process have the opportunity to rescatter. With altered momenta 
in either or both of the decay daughters, reconstruction of the parent particle through 
invariant mass of the daughters is no longer possible. Particles whose decay time scales 
are on the order of the lifetime of the rescattering period (perhaps a few fm/c) might 
then have a suppression in their reconstructed yield proportional to that lifetime. One 
must, however, make some assumption as to the expected yield in order to determine the 
fraction lost. There also remains an argument that some additional production of the 
parent resonances might occur during rescattering. This would dilute the measurement 
of the fraction lost. 

To gain some further understanding, STAR has measured several short-lived species, 
and has more to come. As an example, Fig. | shows the h + tc~ invariant mass spectra 
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Figure 9. The ir + ir invariant mass distributions after background subtraction for the 40%-80% of the 
hadronic Au+Au cross section (left) and for p+p interactions (right) at ,/s NN = 200 GeV. 



from Au+Au and p+p data at y/s NN = 200 GeV. Here we see evidence for such short-lived 
resonances as the i"T*°(892), p(770), and fo(980). Work is also in progress to measure the 
K* +, ~ , A, A(1520), and E(1385) (for the first time in heavy ion collisions), along with 
the 0, which was presented at this conference [ |5|] . 

In particular, the yield of the K*° has been obtained [ |6|] and was already shown in 
Fig. H to not differ notably from a thermal model for central Au+Au collisions. A more 
detailed comparison can be made by taking the ratio of the K*° yield to that of the K 
and c/>. These ratios are shown in Fig. along with the same ratios for other colliding 
systems and energies. On this level, it appears that K*°/K is to some degree suppressed 
in heavy ion collisions at RHIC energies, while the <f)/K*° ratio appears high. Together 
these hint at a slightly low value for the yield of the K*°. This could be construed as loss 
due to in-medium decay. However, with a cr shorter than the size of a Au nucleus (and 
average velocities considerably slower than c), such losses are actually rather small. This 
is consistent with either a very short lifetime of the rescattering period [ |32j , or a loss rate 
only slightly higher than the regeneration rate during rescattering. 

This is only one in a series of observables 
which may help us understand the time 
scales involved in the collision evolution. 
STAR is actively working on measurements 
via other techniques such as HBT corre- 
lations and balance functions to provide 
more handles on time and size scales. Sev- 
eral results from these studies have also 
been presented at this conference [ [33)1 . 




7. Photonic n° Decay 

A notable demonstration of STAR's ca- 
pabilities comes in measuring the spectra 
of 7T° through the 27 decay channel, via 
the subsequent conversion of each 7 into 
an e + e~ pair in the structural material of 

the detectors. Details of the analysis for this 4-body final state can be found in Ref. 



Figure 10: K*°/K and <f>/K*° yield ratios as a 
function of colliding system and energy. For more 
details, see Ref. [ pi. 
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Spectra for the ir° are measured for several centralities. As seen in Fig. [U], the shape 
compares well with the soft negative and positive hadron spectra for central Au+Au col- 
lisions at 130 GeV (which are dominated by 7r~ and 7r + respectively at low pt)- An 
overall normalization uncertainty of ±49% not indicated in the shown errors remains. It 
will be exciting to see what can be done with the improvements available for making 
this measurement in the 200 GeV data and beyond 
from STAR. One addressable topic may be deter- 
mination of the 7T°^77 contributions to the pho- 
ton spectra. This and the 77 — ^77 decay alone ac- 
count for >97% of the inclusive photon spectra in 
Pb+Pb collisions at = 17 GeV [ |4|. STAR 

measures the inclusive 7 spectra (currently via the 
method outlined above using the TPC, and ahead 
with the electromagnetic calorimeter), and has ob- 
served an 77° signal [ [7|; the future looks intriguing. z r«, 
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Figure 11: Spectra of ir° and h + ' from 
central Au+Au collisions at a/snn = 130 



GeV in STAR. Errors shown on 



data 



are statistical and do not include an over- 
all normalization systematic uncertainty of 
±49%. For 7T° spectra from other centrali- 
ties and further details, see Ref. [ pf. 



Data from the STAR Experiment in the soft 
regime at midrapidity probe bulk properties of the 
evolution of heavy ion collisions at RHIC energies. 
Net baryons from these collisions continue to de- 
crease with increasing energy, but are still non- 
zero. This means that pair production processes 
behind baryon yields are now dominant, but do 
not account for all baryon number; there must re- 
main some transport of baryon number even over 
~5 and ~6 units of rapidity at 130 and 200 GeV 
respectively. Also, the (px) systematics of inclusive 
hadron yields provide additional information about 

initial conditions for particle production. Further investigations into photon spectra may 
help us understand the early environment of the collision system. 

Particle yields seem, for the most part, to be well-described phenomenologically by 
statistical models. The data support thermal, chemical equilibrium fits to the yields. Such 
fits proffer a system with thermodynamic properties of temperature and baryo-chemical 
potential in the region of the phase change predicted by lattice QCD calculations. 

Evidence for transverse radial flow is observed in the transverse momentum spectra of 
observed particles. This is consistent with the observation of strong elliptic flow by STAR [ 
35]. In addition, apparently slightly suppressed K*° yields hint at hadronic rescattering 
which is short in duration if regeneration of the K*° resonance is much less than the losses 
caused by rescattering of in-medium decay daughters. 

The STAR Experiment at RHIC continues to address a wide variety of physics topics 
in heavy ion and nucleon collisions. Topics explored through measurements in the soft 
physics regime are strengths of the experiment and will persist in providing a better 
understanding of the collision dynamics throughout the collision evolution. 
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Table 1 : Preliminary results on midrapidity measurements from ^/s NN — 200 GeV Au+Au collisions in STAR. The most central and most peripheral 
measurements available are shown. Inverse slope T, (pr), and x 2 /dof are from mr exponential fits, except for Bose-Einstein on ir~ and pr Gaussian 
on p (| where T is replaced by the Gaussian width a). First error quoted is statistical, second is systematic; if only one error is quoted, the error is 
dominated by systematics. x 2 /dof is large for 7r~, K~ , and p due to point-to-point systematics which are much larger than the statistical errors. 



